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Abstract. By applying the theory of slowly rotating stars to the Sun, the solar 
quadrupole and octopole moments J2 and J4 were computed using a solar model 
obtained from CESAM stellar evolution code (Morel (1997)) combined with a recent 
model of solar differential rotation deduced from helioseismology (Corbard et al. 
(2002)). This model takes into account a near-surface radial gradient of rotation 
which was inferred and quantified from MDI f-modes observations by Corbard and 
Thompson (2002). The effect of this observational near-surface gradient on the 
theoretical values of the surface parameters J2 , J4 is investigated. The results show 
that the octopole moment J4 is much more sensitive than the quadrupole moment 
J2 to the subsurface radial gradient of rotation. 



1. Introduction 

Several theoretical determinations of the J2 and J4 gravitational mo- 
ments have been undertaken in case of different solar differential rota- 
tion laws : (i) only radius dependence (Goldreich and Schubert, 1968; 
Paterno, Sofia, and Di Mauro, 1996), (ii) quadratic expansion in colati- 
tude cosine terms (Ulrich and Hawkins (1981a and 1981b)), (iii) angular 
velocity distribution with a slowly latitude variation determined by 
mean of helioseimology technics (Gough (1982)). More recent determi- 
nations are those performed by : (i) Armstrong and Kuhn (1999) using 
a quadratic expansion rotation law with coefficients obtained by fitting 
higher resolution helioseismic interior rotation data from MDI (Scherrer 
et al. (1995)), (ii) Godier and Rozelot (1999) and (iii) Roxburgh (2001) 
using the differential rotation model given by Kosovichev (1996) from 
BBSO p-modes observations. This model takes into account the pres- 
ence of a constant near-surface radial gradient based on the assumption 
that the angular momentum is preserved in the supergranulation layer. 
The aim of the present work is a contribution to J2 and J4 determina- 
tions using a new analytical model of solar differential rotation provided 
by Corbard et al. (2002) which has a latitudinal dependent profile of 
the near-surface radial gradient of rotation. 

If we consider the Sun as an axial symmetry distribution of matter in 
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rotation, the outer gravitational field 4>out is expressed as : 



4>out{r,0) 



oo , 7-, s 2n 



n=l 



(1) 



where J2 n are the gravitational moments, P 2n the Legendre polynomials 
and r, 6* respectively the distance from the Sun centre and the angle to 
the symmetry axis (colatitude). 

Since the solar rotation is slow, it induces small perturbations around 
the spherical equilibrium. These perturbations can be expanded on Leg- 
endre polynomial basis. The distribution of the gravitational potential 
in the Sun can be written : 

oo 

<fi(r, 6) = Mr) + Mr, 0) = Mr) + E <Wr)P 2 „(u) (2) 

n=l 

where u = cosO. 

The gravitational moments J2 n are obtained assuming the continuity 
of the gravitational potential at the surface : 

R 

The perturbed potential is obtained by linearization of the equation of 
hydrostatic equilibrium and the Poisson equation, leading to : 

d2<t>12n + -^-(2n(2n + l) + UV)^ = U[(V + 2)B 2n + 



dr 2 r dr r 2 

dB 2n 4n + l f+ l 



+ r-^ + ™±± J i (1 - u*)P 2n (u)n(r,u) 2 du} (4) 

where U = A7rpor 3 /M r , V = dlnpo/dlnr, M r is the mass contained in 
a sphere of radius r and Q(r,u) the angular velocity. B 2n is given by: 



B2n{r) = f'u^uf^iiu 2 - !)> (5) 



Equation (3) is integrated with the usual boundary conditions, us- 
ing U and V provided by a solar model obtained from the CESAM 
stellar evolution code (Morel (1997)) and a rotation law derived from 
helioseismology. 
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2. Analytical model of solar differential rotation 

We consider the recent Corbard et al. (2002) model of solar differen- 
tial rotation and, for comparison, the Kosovichev one (1996) already 
used by Roxburgh (2001) and Godier and Rozelot (1999). The Corbard 
model contains a near-surface radial gradient of rotation inferred from 
the radial dependence of the MDI f-modes observations (Corbard and 
Thompson (2002)). Two estimations of this gradient have been derived 
from different sets of modes leading to two rotation models denoted af- 
terward by (a) and (b). As for Kosovichev's model, the surface rotation 
is forced to surface plasma observations (Snodgrass, 1992). 
Figure 1 shows the solar rotation profiles corresponding to these models 
computed for different latitudes. Kosovichev's model presents a nega- 
tive constant subsurface radial gradient. The Corbard models have a 
negative value of the radial gradient at low latitude which is twice 
smaller than Kosovichev ones. At high latitude, the gradients are pos- 
itive with larger magnitudes for the Corbard model (b). 
We use an analytical expression derived respectively by Corbard et al. 
(2002) and Dikpati et al. (2002) for the Corbard rotation laws and for 
the Kosovichev one. We recall hereafter the full set of equations and 
parameters they give for these rotation laws : 

Q(r, u) = Ai (r, u) + ^tac(r) (ft cz - + a>2U 2 + a 4 u 4 ) (6) 

where 

Ai(r,«) = n + y cz (r){a(u)(r-r cz )} + 

+ * s (r){n eq -n cz -f3(u)(r-R Q )-a(u)(r-r cz )} (7) 




Figure 1. Profiles of the solar rotation from O.55i?0 to the surface for different lati- 
tudes computed each 10° from the Equator (top) to the Pole (bottom). (1) Model 
of Kosovichev. (2) Model of Corbard (a). (3) Model of Corbard (b). 
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with 

Q eq - Q cz + f3(u)(R & -r s ) 
a(u) = — 

f2o> ^eq an d ^cz are respectively the constant rotation of the radiative 
interior zone, the equatorial rate at the surface and at the top r cz of 
the tachocline localized at r tac . The and 04 constants describe the 
latitudinal differential rotation. (5{u) represents the latitudinal depen- 
dance of the rotation radial gradient below the surface down to a radius 
r s . This gradient depends on the latitude through the (5q, fa and (3q 
constants by the Equation : (3{u) = (3q + + P&u 6 . 

The VE^ function where x stands for tac, cz or s, models the transition 
between different gradients. An error function centered at r x with width 
lj x is used for this goal : *S> x (r) = 0.5 (1 + erf [2(r — r x )/to x \). 
All the rotation laws have the following common parameters : Qq = 
435nHz, £l eq = 452.5nHz, VL CZ = 453.5nHz, r tac = 0.69i? Q , r cz = 
O.71i?0, 02 = — 61nHz, 04 = —73.5nHz. The parameters which are 
different for the three laws are given in Table I. 



Table I. The non common parameter values between the rotation models. The /3o, P3 
and /3e parameters are given in hHz/Rq. 



model 


OJtac/RQ 


^cz/ Rq 


Us/Rq 


Ts/Rq 


A) 


ft 


ft 


Kosovichev (1996) 


0.1 








0.983 


891.5 








Corbard (a) (2002) 


0.05 


0.05 


0.05 


0.97 


437 


-214 


-503 


Corbard (b) (2002) 


0.05 


0.05 


0.05 


0.97 


437 





-1445 



3. Results and discussion 

We present in Table II the computed values of J2 and J4 obtained 
with different solar rotation models described in Section 2 and with 
an uniform rotation equal to the rotation rate of the solar radiative 
zone CIq, for comparison. Table III gives also some values of J2 and J4 
presented by other authors. 

Our results show that the differential rotation in the convective zone 
reduces J2 value of about 0.8% when the Corbard models are considered 
and about 0.5% in the case of Kosovichev's model. For this last case, our 
J2 determination is larger than the value found by Godier and Rozelot 
(1999) but in agreement with the one obtained by Roxburgh (2001), 
both using Kosovichev's model. Our values are also in agreement with 
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Table II. The J2 and J4 values corresponding to the 
different rotation models. f2 = 2.733 fxrd/s is the 
rotation in the radiative zone. 

Model of rotation J 2 (xlCn 7 ) J 4 (xlCn 9 ) 



Uniform rotation (fio) 
Kosovichcv (1996) 
Corbard (a) (2002) 
Corbard (b) (2002) 



2.217 

2.205 -4.455 

2.201 -5.601 

2.198 -4.805 



those given by Paterno et al. (1996), Pijpers (1998) and Armstrong and 
Kuhn (1999) (see Table III). All these values deviate from the range 
obtained by Ulrich and Hawkins (1981a and 1981b) in the case of the 
rotation law defined as a simple quadratic expansion. The difference 
between the subsurface radial gradients induces only a small reduction 
on J2 values. It is about 0.25% between Kosovichev's model and Cor- 
bard's ones. This difference is however about 0.1% between the two 
Corbard models. 

As expected, the effect of the subsurface radial gradient is more im- 
portant on J4 gravitational moment. J4 absolute values obtained using 
the models of Corbard (a) and (b) are respectively about 20% and 7% 
larger than the one obtained with Kosovichev's model. The (a) Corbard 
model increases the J4 absolute value of about 14% compared to the 
one obtained from the (b) Corbard model. Our \ J±\ value corresponding 
to Kosovichev's model is in agreement with the one given by Roxburgh 
(2001) using the same model. However, those obtained from Corbard's 
models are larger than other values given in Table III. All these values 



Table III. Some computed values of J2 and J4 obtained by other 
authors. The large value of Gough (1982) is due to an estimation of 
the internal rotation deduced from earlier helioseismic observations 

Authors J 2 (xl0- 7 ) J 4 (xl0" 9 ) 

Ulrich and Hawkins (1981) 1.0 < Ji < 1.5 2.0 < |J 4 | < 5.0 

Gough (1982) 36 

Paterno et al. (1996) 2.22 

Pijpers (1998) 2.18 

Godier and Rozelot (1999) 1.6 

Armstrong and Kuhn (1999) 2.22 -3.84 

Roxburgh (2001) 2.206 -4.45 
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are consistent with the range given by Ulrich and Hawkins (1981a and 
1981b) except for the (a) Corbard model. 

Rotation induces a distortion of the solar surface which can be roughly 
related to J2 through the following quantity often called oblateness : 

R e ~ R P _ 3 T , ft 2 s R Q /m 

— R^~2 J2+ 2GMe (8) 

where £l s is an effective rotation rate. R e , R p and R & are respectively 
the equatorial, polar and mean solar radius. This formula is strictly 
valid for an uniform rotation or for a rotation constant on cylinders. 
For a solar rotation which presents a complex profile not constant on 
cylinders, Paterno, Sofia, and Di Mauro (1996) proposed an expression 
of Q s derived from the surface rotation Q(R Q ,u). In our case, Q s will 
be the same for the three models since they are built so as to have 
the same surface rotation. Thus, in this rough description, the effect 
of different subsurface radial gradients of rotation on the oblateness 
appears through the modification of the J 2 gravitational moment. It 
is negligible since the main term in Equation (8) is the surface rota- 
tion term. The value of the oblateness found is 9.1 x 10~ 6 . It is in 
agreement with observations of Lydon and Sofia (1996) and Rozelot 
and Rosch (1997) but slightly larger than those of Kuhn et al. (1998) 
and Armstrong and Kuhn (1999). Rozelot, Godier, and Lefebvre (2001) 
presented new developments taking into account the surface latitudinal 
differential rotation to link J2 and J4 to the solar equatorial and polar 
radius. To the lowest order, their formula reduces to Equation (8) with 
another definition of £l s . For our surface rotation, their formula leads 
to an oblateness value equal to 11.4 x 10~ 6 . We have estimated that the 
effects of the second terms are of the order of (4/1000) in the case of 
Kosovichev's model. New constraints on the oblateness and the shape 
of the solar surface will hopefully be provided by the future PICARD 
microsatellite CNES-mission (Thuillier et al. (2003)). 
In conclusion, the octopole moment J4 is much more sensitive than 
the quadrupole moment J2 to the inclusion of the latitudinal and ra- 
dial differential rotation in the convective zone and particularly to the 
subsurface radial gradient of rotation. Indeed, an important subsurface 
radial gradient at high latitude decreases significantly the value of \J±\ 
while it does not affect significantly the J 2 value. 
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